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Abstract. The pionic contribution to the 
g — 2 of the muon involves a certain in- 
tegral over the the modulus squared of 
F n (t), the pion electromagnetic form fac- 
tor. We extend techniques that use cut- 
plane analyticity properties of F n (t) in 
order to account for present day estimates 
of the pionic contribution and experimen- 
tal information at a finite number of points 
in the space-like region. Using data from 
several experiments over a large kinematic 
range for \t\, we find bounds on the ex- 
pansion coefficients of F n (t), sub-leading 
to the charge radius. The value of one 
of these coefficients in chiral perturba- 
tion theory respects these bounds. Fur- 
thermore, we present a sensitivity analy- 
sis to the inputs. A brief comparison with 
results in the literature that use observ- 
ables other than the g — 2 and timelike 
data is presented. 

1 Introduction 

General principles such as analyticity and unitarity have 
been used to derive useful constraints on form factors (for 
a review see ref. pQ). In particular, useful lower bounds 
were obtained on the pionic contribution to the (g — 2) of 
the muon (muon anomaly). For recent reviews on the cur- 
rent status of this quantity, see refs. [U[3]. Consider the 
following low-energy expansion for the pion electromag- 
netic form factor F n (t), 

F n (t) = F 7T (0) + ^t + ct 2 +dt 3 + --- , (1) 
6 

where rv is the pion charge radius, c and d are Taylor co- 
efficients, with units of GeV~ and GeV -6 respectively. 
For a recent discussion, refer [7J. From general princi- 
ples, discussed in the next section, and using the nor- 



malization i^(0) = 1, immediately gives a lower bound 
of ~ 1.6 x 10 -9 . Imposing the experimental value of r n 
further improves the bound to 2.9 x 10~ 9 [4]. Caprini [5] 
has more recently shown that if one were to assume that 
the pionic contribution is less than I = 75 x 10~ 9 , then 
these principles yield constraints on the higher expansion 
coefficients, c, d which translate into an allowed ellipse in 
the c — d plane. These results, as well as the results we 
present here, hold as long as the true pionic contribution 
is less than the value of / we have assumed. If the pionic 
contribution is indeed significantly smaller, then the true 
allowed region would be a proper subset of the region we 
isolate. 

The inclusion of experimental data from the time-like 
region, the phase, the modulus, and the phase as well 
as the modulus each improves the bound on the muon 
anomaly pQ. Using / as an input, Caprini has used dis- 
persive techniques and time-like data, and has considered 
other physical observables such as certain QCD polariza- 
tion functions, in addition to the muon anomaly, to con- 
strain further the allowed regions in the c — d plane Q. We 
note that the parameter c itself has been determined to 
2-loop accuracy in chiral perturbation theory [7] and has 
a value of 4.49 GeV -4 , which is safely accommodated in 
the allowed ellipse of Caprini [5]. 

An important source of experimental information are 
the values of the electromagnetic form factor in the space- 
like region. These have been reported in electroproduction 
experiments conducted in the 1970's (see refs. [5], [9], [TU]). 
Raina and Singh [llj used this information to produce im- 
proved bounds on the muon anomaly (see also |12j). New 
experimental data is now available (see refs. [T3], [II], [T5] 
and [16]). In this work, we will make use of some of the 
data sets above. 

Our main purpose in this work is to demonstrate that 
the framework of Raina and Singh can be effectively used 
to obtain constraints in the c — d plane, using the pionic 
contribution to the muon anomaly as an input. This is an 
algebraic framework, which is clear cut and transparent. 
It provides an important consistency check on the allowed 
regions isolated by Caprini. 

In sec. [2] we briefly review the dispersive formalism 
and describe the implementation of space-like constraints 
and the method by which we isolate the region in the c — d 
plane. Here, we provide the theoretical framework and also 
a discussion on the present experimental status, as they 
are both required for this purpose. In sec. [3] we present 
our results using data from [8], [9], [10] and more recent 
data from [T5] and [IB]. We present a detailed discussion 
of our results and our conclusions in sec. [4] 



1 Applications to form factors in semi-leptonic decays of 
such general methods in a modern context have also been re- 
cently considered [6]. 
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2 Formalism 

We recall the formalism presented in ref. [UJ. The pion 
contribution to the muon anomaly is given by: 



dt p(t)\F w (t)f 



(2) 



where t„ = Am 2 , is the branch point of the pion form factor 
and 

a 2 m 2 (t - Ufl 2 
P(t)= J 1 t yl K(t)>0 (3) 



12tt 



where, 



K{t) 



2 2 
771 U 

du(l -u)u 2 {\ -u+ — ^— (4) 



Using the following map from the i-plane, which is cut 
from along the real t axis, to the complex z-plane (re- 
gion \z\ < 1), 

z 1 . jt t-Tj- 

(5) 



2 + 1 



and the definitions: 



/(*) = F*(t) 
p(z) = p{t) 



(6) 
(7) 



the pionic contribution to the muon anomaly can be writ- 
ten as 



1 



2tv 



a,(n+n-) = ^- J dO w(e)\f(e iS )\ 2 
where, 

w{0) = Am\ scc 2 (6>/2) tan(0/2) p (e ie ) > 0. 
We now consider a function h(z) defined as: 

where, 



w n (z) = exp 



1 f 2 * e" 
-— / d0 6 

47T 



Then eqn. 1(5)) can be written as: 



a^ir+n-) = — / M\h(e") 
27T Jo 



i9\\2 



(8) 

(9) 
(10) 

(11) 
(12) 



Now h{z) is analytic within the unit circle \z\ < 1 and 
for real z, h(z) is real. Therefore h(z) can be expanded as 
follows: 

h{z) = ao + aiz + ci2z 2 + • ■ • (13) 

where ao, &i • • • are real coefficients. Therefore, in the an- 
alytic region, \z\ < 1, a M (7r + 7r~) can be written as: 



a M (7r + 7r ) = al + a\ + 



(14) 



The expansion coefficients a n can be obtained from a Tay- 
lor expansion of the function h(z) in terms of f(z) and 
w 7r (z). The coefficients are given by 



a = h(0) = w„(0), 
2 



(15) 



a a = fc'(O) = <(0) + ^t w ttf w (0), (16) 



(l 2 



h"{0) 



2! 



^(0) 
( \rlt« 



32 ctl 



<(0) 



(17) 



and 



/i"'C0) 1 

a 3 = -^p = g M°) (l2t*U - 384 ct* + 384 dtl)} 
1 

+ 6 



3<(0) ( --r*t w + 32c# 



+ -[2<(0)r^ + <"(0)]. 
In our treatment, the expansion coefficients satisfy 



i=0 



(18) 



(19) 



Given /, eqn. (|19p yields constraints on the expansion co- 
efficients of the form factor. Including up to the second 
(third) derivative for F w (t) results in constraints for c (c 
and d). It may be pointed out that Caprini's result on c 
and d shown in Fig. 1 (dashed lines) of ref. [3] is obtained 
precisely in this manner. 

The constraints on the expansion coefficients of inter- 
est may be significantly improved through the inclusion of 
experimental information on the form factor. Our objec- 
tive in this work is to study the effect of including experi- 
mental information coming from the space-like region. In 
order to meet this objective, we first extend the formalism 
that has been presented in ref. [TT]. In that work, space- 
like constraints were used to obtain lower bounds on the 
muon anomaly, denoted as I m i n . 

We observe that N space-like constraints are linear 
constraints that may be expressed as: 



a,c 



(») 



(20) 



i=0 



where n = 1, 2, ■ ■ • N. Such constraints can be implemented 
through the technique of Lagrange multipliers, by setting 
up the Lagrangian 

oc N oo 

i = oE a ' + E^«-E ) ( 21 ) 



i=0 n=X 



i=0 



where in eqn. (|2ip . we consider only finite number of ex- 
pansion coefficients. The Lagrange equations yield: 



E 



a n c. 



(n) 



(22) 
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and 



N oo 



» 



It then follows that 



n = l i=0 



N 



n=l 



where 



to these experimental uncertainties as their dimensions 
increase. This sensitivity is particularly severe for small 

(23) 

values of \t\. As a result, at such values we are able to 
implement at most two constraints, while data from higher 
energies allows us to implement up to three constraints. 
However, the data from smaller values of \t\ provide more 
stringent bounds for fixed number of constraints. Thus we 

(24) 

see a fairly complex interplay between the energy regime 
that we can use and the number of constraints we are able 
to implement. 



(25) 



N oo 
m— 1 i— 

All the a„s may be eliminated to yield a determinantal 2.1 Space-like data 

equation for / 



I A x A 2 



(26) 



For the case at hand, where we wish to specify ao, ai , a 2 , CI3 
and the value at space- like points h{x,i), i = 1, 2, 3..., where 
Xi is real, the determinantal equation reads, 





a ai a 2 


«:-! 


h(xi) 




h(x 2 ) 


"(I 


1 





1 
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a x 


1 









x 2 
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1 
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03 
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x\ 




x 2 


h(xi) 


1 Xi x\ 


x\ 


(1-x?)- 1 


(1- 


- XiX 2 


h(x 2 ) 


1 x 2 x\ 


x\ 


{l-X 2 Xi) 1 


(1 


-A) 



(27) 

In the above merely retaining the first two rows and columns 
gives the first bound of Palmer, given in the Introduction, 
while retaining the first three rows and columns yields 
the second of Palmer's bounds. Raina and Singh [11] use 
the value of ^(0) = 1 and rv to obtain the lower bound 
to a fl ('!T + TT~). This amounts to dropping out the rows and 
columns corresponding to a 2 and 03 in the eqn. ([27]) which 
are related to the expansion coefficients c and d. Instead, 
providing an input to / in eqn. (|27p . gives us an allowed 
region in the c— d plane. Dropping the row and column cor- 
responding to 03 would result in determining an allowed 
region for c alone, which we pursue in the next section for 
purposes of illustration. 

In the next subsection, we provide a discussion on the 
present day experimental information that is utilized in 
our study. This information spans an impressive range of 
energies, viz. \t\. While in principle there is no limit to 
the number of constraints, in practice the uncertainties 
in the experimental determination and sensitivity of the 
determinantal equation limits this number. The reason for 
such restrictions is the extreme sensitivity of the matrices 



We will begin with the data in the space-like region that 
was used in the work of Raina and Singh [11] , that came 
from measurements in the seventies. We shall refer to these 
as the Brown data and Bebek data respectively. The data 
we use from these sets is given in Tables [1] [2j The ta- 
bles also list the values of z(t) = x(t), which is the map 
from the t plane to the disc \z\ < 1 (x(t) lies in the range 
[—1 < x(t) < 0]), for the chosen data points and the corre- 
sponding value of h(z) = h(x) (refer eqn. (fT0|) ). To make 
our notations clear, we refer to the data points correspond- 
ing to a particular \t\ as x\, x 2 and so on, in the ascending 
order of magnitude of \t\. The tables also show the exper- 
imental errors in the data. 

More recent data come from four sources. The first 
of these is ref. |13j which is from a Fermilab experiment 
(Fl). Here the range is 0.03 GeV 2 < -t < 0.07 GeV 2 . 
Even at the high end, we are way below the lowest en- 
ergy of Brown. The next set of data come from another 
Fermi lab experiment (F2) (ref. [H]), covering the range 
of 0.037 GeV 2 < -t < 0.094 GeV 2 . The CERN NA7 ex- 
periment provides very accurate data [15] in the range 
0.014 GeV 2 < -t < 0.26 GeV 2 (by Amendoliaet.al). There 
is an overlap region between this last set and the Brown 
data and the two are consistent. The final data we use 
come from the JLab experiments (Tadevosyan et ah, ref. [16]). 

In our study we choose to work with recent data from 
Amendolia et.al [T5] that cover both the low and interme- 
diate energy range (note that the data from [T3lfT4] cover 
a smaller energy range; hence our choice) and the data 
from the JLab experiments (Tadevosyan et al., ref. [H]) 
that cover a higher range of \t\. Using really small \t\ val- 
ues lead to numerical instabilities as the entries in the 
determinant (eqn. ([27]) ) become small. 

Tables [T] [2] [3] [4] show the data that is used in our 
analysis to constrain the expansion coefficients c and d 
of the form factor F w (t). The data from various experi- 
ments are chosen so as to give reliable numerical results. 
We find data in lower \t\ to constrain the expansion coeffi- 
cients better but we cannot include more than two space- 
like constraints; while the higher \t\ region gives a weaker 
bound, which can be improved by increasing the number 
of space-like constraints. 
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Table 1. Space-like data from Bebek et.al [TO] 



f(-Q 2 ) [GeV 2 ] F n (t) x(t) h(x) x 10" 5 

1 -0.620 0.453 ± 0.014 -0.499 3.057 

2 -1.216 0.292 ± 0.026 -0.606 2.035 

3 -1.712 0.246 ± 0.017 -0.655 1.716 



Table 2. Space-like data from Brown et.al [9] 



t(-Q 2 ) [GeV 2 ] F v (t) x(t) h(x)xW~ 5 

1 -0.294 0.606 ± 0.028 -0.372 3.775 

2 -0.795 0.380 ± 0.013 -0.540 2.608 



Table 3. Space-like data from Tadevosyan et al. [16] 



t(-Q 2 ) [GeV 2 ] F n (t) x(t) h(x) x HT 5 

1 -0.600 0.433 ± 0.017 -0.494 2.915 

2 -1.000 0.312 ± 0.016 -0.576 2.163 

3 -1.600 0.233 ± 0.014 -0.645 1.626 



Table 4. Space-like data from Amendolia et.al [15] 



i(-Q 2 ) [GeV 2 ] 


F w (t) 


x{t) 


h(x) x 10~ 5 


1 -0.131 


0.807 ± 0.015 


-0.242 


4.454 


2 -0.163 


0.750 ± 0.016 


-0.275 


4.286 



3 Results 



In our results, we use I = 75 x 10 -9 as an input, as quoted 
in [5]. This value is certainly greater than a recent esti- 
mate for the hadronic contribution to a!i ad = 69.2 x 10~ 9 
given in [3] . Pionic contribution is expected to be around 
70% of the total hadronic value. Our value for / gives a 
conservative upper bound to the pionic contribution to a M . 

We begin by constraining only the expansion coeffi- 
cient c which amounts to dropping out the row and col- 
umn corresponding to 03 in eqn. ([27]) . We determine the 
bounds on c in stages, starting with only the pion charge 
radius r w and the normalization of F„(t) at t = (i.e., 
no space-like constraints) and then incorporate space-like 
constraints successively. The results arc shown in fig. ([T]), 
where we have used data from tables [T] [2] [3] and [4] For 
any data set, we start with the smallest value of | £ | for 
one space-like constraint and include more constraints in 
increasing magnitude of \t\. All the data sets show the 
trend that inclusion of more space-like constraints in this 
manner, improve the bounds on c. The largest range for 
c in units of GeV~ 4 is obtained when no space-like con- 
straints are used. This agrees with the result of Caprini 
(fig.(l) in [5]). Note however, that for data from larger 
values of \t\, the bound on c is significantly weaker, com- 



pared to data from smaller values of \t\ for a fixed number 
of space-like constraints. 

It is worth investigating the sensitivity of the bounds to 
the errors in the data. The Brown data shifts the bounds 
on c to negative values when two space-like constraints are 
used. This shift is sensitive to the error bounds on F^(t). 
If we use as input F v {t x ) + 0.028 (table [2]) for the first 
constraint and F^fo) — 0.013 for the second constraint, 
then the upper bound on c is positive. Here, t\ and ti 
represent the first and second data point respectively in 
table[2] We present further sensitivity analysis later in this 
section. Similarly for the data from Amendolia et.al |15| . 
we see that one space-like constraint can be incorporated 
using central values for F n (t) ; but varying this input over 
the error bounds given in table [4] we can use up to two 
constraints (not shown here). Unless othewise specified, 
all results presented in this section use the central values 
of the form factor data given in tables Q] [2] [3] and [4]. 

Including constraints from 03 results in a relationship 
between the expansion coefficients c and d which we shall 
now explore. Fig. ([2]), shows the allowed region in the c — d 
plane using only the normalization condition on F^{t) and 
the value of rv . Analogous to the study of the bounds on c, 
we can include space-like constraints to impose stringent 
bounds on the expansion coefficients c and d, that are 
reflected by smaller allowed regions in the c — d plane. 
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Fig. 1. Constraints on c starting with just information on F^(0) = 1 and r n and incorporating space-like constraints. Note that 
the allowed values of c decrease as we include more space-like constraints. 
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Fig. 2. Constraints on c and d using the normalization F^{Q) = 1 and the pion charge radius 7V. 
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value of h(z) — h(xx) and study the variation of the bounds on c and d as we use data from higher the \t\ region. 




-40 -20 20 -30-20-10 10 -10 10 



> 

CD 

o 




-40 -20 20 
c [GeV" 4 ] 



30-20-10 10 
c [GeV" 4 ] 



c [GeV" ] 



Fig. 4. Constraints on c and d as the number of space-like constraints are increased. Here the space-like data is taken from 
refs. [10], [16]. 
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Fig. 5. Constraints on c and d as the number of space-like constraints are increased. Here the space-like data is taken from 
refs. [9], [H]. 



Fig. ([3]) shows the variation in the bounds on c and d 
for one space-like constraint, where data from different \t\ 
region is used. In the figure, h(x\) is the value of h(x) as 
defined in cqn. (|10[) corresponding to the smallest value of 
t in each data set (refer tables [T]- [4j) . Similarly h{x2) is the 
constraint at a higher \t\ value taken in the ascending order 
of magnitude. It may be readily observed that the most 
stringent bounds are obtained when data from smallest 
values of \t\ are used, and they improve as the number of 
constraints are increased, as noted earlier. 

Fig. ((4| (as well as fig. (0) shows that increasing the 
number of space-like constraints improves the bounds on 
the expansion coefficients c and d. We observe that the 
allowed range for c corresponds to that in fig. [1] Using 
the central values tabulated in tables [1] 0] we see that 
we are able to use up to three space-like constraint with 
the data from Bebek et al. [10] and the recent data from 
Tadevosyan et al. [16], while we are able to incorporate 
only two space-like constraints with the data from Brown 
et al. [pj and one with the data from Amendolia et al. [T5] . 
The best estimate for the bounds on the Taylor coefficient 
c and d is obtained for the data from Bebek et al. [TP] . 
But varying the Bebek data within the error bounds, we 
see that there is overlap with the bounds obtained with 
the data from Tadevosyan et al. [16] , as seen in fig. ([6J . In 
fig. ©, the label "max" refers to the data taken from the 
tablcs[T]and[3]with the corresponding error bounds added, 



i.e. at the upper limit of the error bound, central refers to 
the central values and "min" refers to central value minus 
the error bound, i.e., the lower limit of the error bound. 
We see that an overlap between the bounds obtained from 
the different data set occurs when the data are close to the 
upper error bound. In fact at the upper error bound, the 
data from Tadevosyan et al. gives better bounds on the 
expansion coefficients compared to the data from Bebek 
et al. At the central value of F„(i), the data from Bebek 
et al. does better while at the lower end, the data from 
Bebek et al., no longer yields stable result, while the data 
from Tadevosyan et al. still gives reasonable bounds. We 
find that despite small differences in the allowed regions, 
the fact that the allowed regions are essentially the same 
offers an important consistency check on the form factor 
determinations by each of the experiments. 

As seen from the our results so far, the bounds on 
the Taylor coefficients vary due to the errors in the data 
presented in the tables [1] - S] hence it is important to 
examine the sensitivity of our results to the inputs. We 
have made a systematic study of the sensitivity by varying, 
(a) the input value of rv, (b) the input value of /, and (c) 
varying the experimental determinations of F v (t) within 
their quoted errors, and in each instance keeping all other 
inputs fixed. 

Variations in the bounds on c and d for one space-like 
constraint from [10j . [15] , when the value of rv is var- 
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Fig. 6. Variations in the bounds on c and d as data from [10J and [16] are varied within their error bounds for three space-like 
constraints. The shaded region represents the bounds from [16] and the symbol represents those from |10| . 



ied within the error bounds quoted therein, is shown in 
fig. ^j}. We see that the bounds do not vary much as 
TV is varied within the allowed errors of a few percent. 
Fig. ([5]) shows the variations in the bounds of c and d as 
the value of / is changed. Lower / values results in a sub- 
space of the region allowed for the coefficients at higher 
/ values. We expect this behavior as our method gives 
only an allowed region for the coefficients. This explicitly 
demonstrates that as long as the true pionic contribution 
is lower than the values we have used, the allowed region 
lies within the isolated ellipse. 

Fig. ([9]) shows the variations bounds on c and d for 
the data from Amendolia et.al [TS] and Bebek et.al [TU] 
for one space-like constraint (smallest magnitude of |i| in 
the tables [1] and [4]). The data are varied within their ex- 
perimental bounds and the results arc depicted in fig. §§§ . 
We see that the bounds on the expansion coefficients are 
most sensitive to errors in the data from lower \t\. Hav- 
ing said that, it is our view that one may read-off reliable 
ranges for both c and d in light of this sensitivity analysis. 
However an analogous sensitivity test that we have carried 
out for two space- like constraints for the data from Amen- 
dolia et al. [15] leads to a complete loss of coherence. This 
is an unavoidable numerical difficulty in the determinant 
that arises due to the fact that the entries are small and 
closely spaced. Therefore we can conclude that for low \t\ 
information we are unable to obtain reliable results for 
more than one constraint. 



We explore this sensitivity of the bounds to the data 
taken from the low \t\ region using one space-like con- 
straint from [15j and varying this constraint over the given 
error bounds (Tabled and study the corresponding vari- 
ations in the bounds on c alone. This is seen in fig. (|10[) 
where the circles represent the maximum and the squares 
represent the minimum of the bounds. The filled symbols 
stand for the bounds obtained from the data at lower value 
of | i | and the open symbols represent those from the data 
at higher value of \t\. Fig. ([TTj) shows similar variations 
in the bounds on c when one as well as two space-like 
constraints (ref. table [3]) for the data taken from [16] are 
varied over the error bounds. Note that the data from [16] 
lie at higher \t\ value compared to the data from [IB] . 
When two space-like constraints are used, we fix one con- 
straint at the central value and vary the other over the 
error bounds quoted in [16] and plot the maximum and 
minimum of the bounds as a function of the constraint 
that is varied. From figs. ([TO]) and (|11|) . we see that data 
from lower \t\ value in each set constraints the coefficients 
better for one space-like constraint. Also the sensitivity to 
the error bounds is greater for a set at lower |i| compared 
to the set at higher value of \t\. We also note that two 
space-like constraints are more sensitive to the variations 
in F v (t) compared to one space-like constraint, as can be 
seen in fig. [TO 

Relative errors in the bounds on c as one space-like 
constraint is varied over the quoted error bounds com- 
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pared to the bounds obtained from the central value is 
shown in fig. (| 1 2[) for the data sets from Amcndolia et 
al. and Tadevosyan et al.. We see that relative errors are 
smaller for the data from Tadevosyan et al., once again 
emphasizing our earlier observation that the sensitivity of 
the Taylor coefficients to the errors in the determination 
of F v (t) is greater for the data from the low \t\ region 
compared to those from the high \t\ region. 
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bounds. Note that the bounds are more sensitive to the errors 
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4 Discussion and Conclusions 

In this work we have considered the constraints on the pion 
electromagnetic form factor coming from present day es- 
timates of the pionic contributions to the muon anomaly, 
and from space- like data which is available over a fairly ex- 
tended kinematic regime. We have adopted the framework 
of Raina and Singh and have extended it to obtain con- 
straints on c and d, the expansion coefficients sub-leading 
to the charge radius of the pion. We have used data avail- 
able in the 1970's (also used in ref. [11]) and more recent 
data which spans significantly lower values of \t\, as well 
as very recent data coming from Tadevosyan et al. which 
are at higher \t\ values, which lie in the range of Bebek et 
al. 

As mentioned in the Introduction, Caprini [5] has used 
timclike data to obtain constraints on the Taylor coeffi- 
cients. By systematically considering the inclusion of only 
the phase, the ellipse with no constraints was found to 
shrink considerably using as an observable the QCD po- 
larization function (see figs.(l) and (2) in [5j). This anal- 
ysis employs an optimal technique, resulting in an inte- 
gral equation of the Frcdholm type that was solved nu- 
merically. The ellipse was found to shrink even further 
when both the phase as well as the modulus information 
were used even though the technique was non-optimal (see 
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fig. (3) in [5]). Questions remain about its validity due to 
certain mathematical difficulties, as discussed by Caprini. 

The best estimates for the bounds on c and d is ob- 
tained using the data from Bebek et al. [10] for three 
space-like constraints. In fig. (|13|) . we compare our best 
result with those obtained by Caprini using only the phase 
from the timclike region. Using space-like data we get the 
following range for the expansion coefficients: —1 GeV" 4 
< c < 11 GeV" 4 and -132 GeV" 6 < d < 220 GeV" 6 . 
On the other hand Caprini has a range of [—14 GeV" , 
44 GeV" 4 ] and [-236 GeV" 6 , 594 GeV" 6 ] for c and d 
respectively using only the phase of timclike data. It is 
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interesting to note that the overlap region between the 
two determinations, ours using space-like data and the 
muon anomaly and Caprini's using the phase of timelike 
data and the QCD polarization observable, accommodates 
comfortably the value of c from chiral perturbation the- 
ory. We note here that similar bounds obtained using the 
central values of the Tadevosyan data (ref. table [3]) us- 
ing three space-like constraints does not accommodate the 
value of c obtained from chiral perturbation theory, but 
varying the data at the upper end of the error bound, as 
in fig. [6] the allowed region does include this value. 

Our conclusions are that the constraints are weaker 
than those obtained by Caprini using both phase and mod- 
ulus of data from the timelike region (see fig. (3) in [5]), 
while they are more stringent than those obtained with 
only the phase of the timelike data. Here we note again, 
that the latter of these treatments is rigorous and the re- 
sults may be taken as reliable. On the other hand, the 
results obtained from using both the modulus and phase 
while appearing more stringent are on a less rigorous foot- 
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ing, due to inherent mathematical difficulties as noted in 
ref. 0. 

We have also carried out a sensitivity analysis by vary- 
ing the estimate for the pionic contribution over a signif- 
icant range, varying the charge radius over its presently 
known errors, and also by varying the experimental data 
that we use over its errors, for the values of \t\ that we 
have chosen. Our conclusions are that in the small \t\ re- 
gion, the system is very sensitive and as a result, we have 
been unable to implement anything more than one space- 
like constraint. For larger values, we are able to include 
two and even three constraints. 

Of related interest, is the extreme sensitivity of the sys- 
tem when bounds are being derived on the muon anomaly 
from space-like data, as recognized earlier by Pantea and 
Raszillier [T7] . The issue was of some significance because 
several space-like constraints were being used to obtain 
these bounds. Our circumstances are somewhat mitigated 
by the fact that we are using significantly smaller num- 
ber of constraints, as the determinantal equations we are 
solving are already of rather large dimensions. By study- 
ing the sensitivity of the bounds, we are confident that 
our results remain stable. 

It is conceivable that as data improves one may obtain 
better constraints on c and d, where it is numerically fea- 
sible, that can then be used as inputs for high precision, 
self-consistent form factor determinations. We also note 
at this point that our work takes into account only the 
information present in the space-like region. It would be 
worth-while to properly formulate the problem so that the 
information available in the timelike region could be in- 
corporated together with a comprehensive error analysis. 
A theory of error functionals has been developed by Raina 
and Singh [18] in the context of finding a lower bound for 
the muon anomaly, that could in principle be extended 
for the problem at hand, which is beyond the scope of 
this work. Therefore, future investigations could combine 
highly accurate phase and modulus timelike information 
coming from recent experiments and the techniques devel- 
oped by Caprini for the phase problem, together with a 
suitable extension of the error functional method, to pro- 
duce stringent constraints on c and d. 
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